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We  consider  the  nonlinear  Volterra  integrodifferential  equation 
t  3  .  . 

-py  =  /  a(t-s)  (^ -  -  y  (s)  )ds  -  f(t)ya 

-oo  y  (s) 


N  -X.u 

where  a(u)  =  £  K.  e  1  ,  0  <  a  <  3,  and  f:  3R  -*-]R  is  a  given  function. 

i=l 

We  study  the  existence  of  solutions  and  their  asymptotic  behavior  as  t 
for  certain  classes  of  functions  f.  The  main  result  is  that  for  each  f 
satisfying  t  lira ^  e  CTtf(t)  =  0  for  some  a  >  0  and  f(t)  =  0  for  t  S  tQ 
there  exists  a  unique  solution  y(t)  satisfying  lim  y(t)  =  1.  Moreover, 
for  this  solution  lim  y(t)  exists.  This  holds  both  for  V  ?  0  and 
y  =  0.  It  is  further  shown  that  the  condition  f (t)  =0  for  t  >  tQ  can 
not  be  replaced  by  exponential  decrease  unless  f  is  small  or  a  =  1. 
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SIGNIFICANCE  AND  EXPLANATION 

V 

The  equation  describes  the  evolution  of  the  shape  of  a  filament  or  a 
sheet  of  a  polymeric  liquid  subjected  to  an  external  force  f(t).  y  denotes 
the  length  of  the  filament  or  the  thickness  of  the  sheet  respectively,  and 
u  is  a  viscosity  constant  which  can  either  be  positive  or  zero.  The 
exponent  a  depends  on  the  physical  situation  that  is  being  considered. 
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We  consider  the  length  at  t  =  as  known  and  investigate  the  evolution. 
It  is  shown  that  for  a  physically  realistic  class  of  functions  f  there 
exists  a  unique  solution,  and  that  the  length  approaches  a  new  stationary 
value  at  t  =  00  (which  is  in  general  different  from  the  value  at  t  =  -«) . 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


EVOLUTION  OF  THE  SHAPE  OF  A  POLYMER  SUBJECTED  TO  A  FORCE 


M.  Renardy 

1.  INTRODUCTION 

In  [1]  Lodge,  Me  Leod  and  Nohel  have  studied  the  history  value  problem 
for  the  nonlinear  Volterra  integrodifferential  equation 

t 

-yy(t)  »  /  a(t-s)  F(y  (t),y  (s)  )ds 

—  00 

This  equation  arises  as  a  mathematical  model  for  the  stretching  of  a  filament 
or  a  sheet  of  a  molten  polymer  which  is  assumed  to  be  spatially  homogeneous. 
In  this  case  y  denotes  the  length  of  the  polymer,  y  is  a  parameter 
related  to  the  viscosity, 

3 

F(y,z)  =  X-t  -  z 
z 

and 


N  -A.  u 
a  (u)  =  T  K.  e  1 

■  i  i 
i=l 

with  certain  positive  constants  1C  and  A^.  It  was  assumed  in  [1]  that, 
for  t  £  0,  y(t)  =  g(t)  was  given,  g(-“)  =  1,  and  g  nondecreasing.  A 
general  class  of  functions  a  and  F,  which  includes  those  mentioned,  was 
considered.  One  of  the  results  of  [1]  is  that  under  these  conditions  the 
history  value  problem  has  a  unique  solution,  which  is  nonincreasing  for 
t  £  0  ,  and  approaches  a  constant  strictly  greater  than  1  as  t  -*■  ». 

This  paper  deals  with  a  closely  related  problem.  We  consider  the  elongation 
of  a  polymer  satisfying  the  same  constitutive  law,  but  rather  than  pre¬ 
scribing  the  length  history  we  prescribe  the  force  acting  on  the  polymer, 
which  causes  the  elongation,  and  consider  only  the  length  at  t  =  -®  as 
known.  Unlike  [1],  we  explicitly  use  the  form  of  F  and  a  as  given  above. 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041  and 
the  Deutsche  Forschungsgemeinschaft. 


It  can  be  derived  (for  an  explanation  of  the  physical  principles  see  [2]) 
that  the  elongation  of  the  polymer  is  now  described  by  the  equation 


-uy  =  J  a(t-s)  (^2  —  -  -  y(s))ds  -  f(t)ya  (1.1) 

-00  y  (s) 

The  exponent  a  depends  on  the  physical  situation.  One  case  of  interest 
is  that  of  a  filament  pulled  at  its  ends 


f 


y 


f 


In  this  case  a  =  2.  Another  physically  interesting  situation  is  a  sheet 
of  polymer  pulled  in  two  directions. 


-f 


In  this  case  y  denotes  the  thickness  of  the  sheet  and  a  =  j.  Our  analysis 
holds  in  general  for  0  <  a  <  3.  The  main  result  we  are  going  to  prove  is 
that  for  each  f  satisfying  _  lim  e  0tf (t)  =  0  for  some  a  >  0  there 

t  ">■  — 00 

exists  a  unique  solution  of  the  equation  satisfying  .  lim  y(t)  =  1. 

(Theorem  3.1.).  Moreover,  if  f  =  0  for  t  greater  than  some  t  <  <*>, 

0 

then  this  solution  approaches  a  constant  as  t  -*•  +«>.  Counterexamples  show 
that  this  last  condition  on  f  cannot  be  replaced  by  exponential  decrease 
as  t  -*•  +®  except  for  small  f  (Theorem  2.4.). 
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2.  SOLUTIONS  FOR  SMALL  FORCES 

N  -A.u 

Let  us  consider  equation  (1.1),  where  0<a<3,a(u)=  £  K. e  1  , 

i=l 

and  y  >  0.  This  equation  can  be  reduced  to  a  system  of  ODE's  in  the  follow¬ 
ing  way.  We  put 


t  -A.  (t-s)  . 

Ji(t)  "  /  Kie  1  T~  ds 

y  (s) 


t  -A. (t-s) 

hi(t)  =  /  ^e  1  y(s)ds 


Then  (1.1.)  reads 


-yy  =  l  (g.y  -  h  )  -  f(t*,/a 
i=l 

Ki 

*i"  -Xigi+~ 


(2.1.) 


hi  =  'Xihi  +  Kiy 
2  ^i 

If  we  put  y.  *  g.y  ,  6.  =  —  ,  we  obtain 
i  1  i  y 


n 

-yy  =  l  (y,  -  M  y  -  f(t)y 

i=i 


Y.  =  ~ A . y .  +  K,  -  — Y.  T(y.  -  5.)  +  ^Y.f(t)ya  1  (2.2.) 

1  ii  1  y  i  “  j  j  yi 


6.  =  -A. 6.  +  K.  +  —6.  ^ (y .  -  «.)  -  is.f(t)y 
x  ii  x  j  j  yx 


a-1 


Both  forms  (2.1.)  and  (2.2.)  will  be  used  in  the  following.  Clearly,  if 

K. 

f  =  0,  we  have  the  stationary  solution  y  =  1,  g.  =h.  =  —  =  y.  =  6, . 

a  i  i 
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LEMMA  2.1  The  linearization  of  (2.1)  (or  (2.2)  at  the  stationary  solution 

K. 

y  =  1,  g.  =  h.  =  —  has  zero  as  a  simple  eigenvalue.  All  the  other  eigen- 

1  1  i 

values  have  negative  real  parts. 

PROOF;  Clearly,  (2.1)  and  (2.2)  give  the  same  eigenvalues.  Let  us  consider 
(2.1).  We  obtain  the  following  matrix  of  the  linearization 


N  3K. 

y  1 

_  1 

_  ^ 

_  1_ 

1_  1 

1 

ii  ^xi 

M 

u 

V 

u  u 

V 

-2Ki 

-X1 

0  . 

.  .  o 

0  0  . 

.  .  0 

-2k2 

0 

"X2  * 

.  .  0 

0  0  . 

.  .  0 

-2kn 

0 

0  . 

.  .  -X 

n 

0  0  . 

.  .  0 

K1 

0 

0  . 

.  .  0 

-X1  • 

.  .  0 

K2 

0 

0  . 

.  .  0 

0  -x2. 

.  .  0 

kn 

0 

0  . 

.  .  0 

0  0  . 

• 

This  yields  the  characteristic  polynominal. 


3K. 


3K. 

l 


pm  »- 1 

j  J  ii  ii 


Thus  N  eigenvalues  are  given  by  X  =  -X^  ,  the  remaining  N  +  1  eigen¬ 
values  are  the  zeros  of  the  last  factor.  Obviously  one  of  these  is  zero, 
and  it  is  simple.  It  remains  to  be  proved  that  all  the  remaining  solutions 
have  negative  real  parts.  Assume  the  contrary,  i.e.  s  0,  X  i*  0  and 


3K. 

-  J  — 

i  yXi 


3K. 

1 


l  U(-X.-X) 


0 
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When  we  take  the  real  part,  this  implies 


3K. 

-  I  JT  -  Re  1  +  I 
11 


3K.  (X.+ReX) 
1  i 


i  y[(X^+ReX)2  +  (ImX)2] 


=  0 


It  is,  however,  easy  to  see  that  the  left  hand  side  is  strictly  negative, 
whence  we  have  a  contradiction. 

Let  cj  denote  the  space  of  all  T-periodic  Cn- functions  IR  -*•  EL  The 
following  lemma  holds: 

LEMMA  2.2  If  f o  (t)  £  c£  has  sufficiently  small  norm,  then  there  exists  a 

one-parameter  of  constants  C  and  functions  Y  =  (y,Yi»|5^)  £  (C^+^)2N+1 

such  that  Y  satisfies  equation  (2.2)  with  f  *  f  -  C. 

PROOF:  Obviously  (2.2)  with  f  =  f  -  C  can  be  written  in  the  form 

G(Y,c,f  )  =  0,  where  G  is  a  smooth  function  mapping  a  neighborhood  of 
U  K  K 

„  ,,  i  i.  .  .  ,„n+1.2N+l  _  n  .  .  ,_n. 2N+1 

Y  =  Yq  =  d,  — ,  — ) ,  c  =  0,  f0  =  0  in  (CT  )  xm  x  CT  into  (CT) 

It  is  a  straightforward  consequence  of  lemma  2.1  that  the  Frechet  derivative 
D^G (Yg,0,0)  has  a  one-dimensional  kernel  and  its  range  has  codimension  1. 
Moreover,  one  easily  sees  from  (2.2)  that  DcG(YQ,0,0)  is  not  in  the  range 
of  D  G(Y  ,0,0),  i.e.  D.  .G(Y  ,0,0)  is  onto  and  has  a  one-dimensional 

I  U  V*  t  C  /  U 

kernel.  The  lemma  now  follows  from  the  implicit  function  theorem. 

DEFINITION  2.3  Let  X°  =  { f  e  OR,  3R)  |  lim  e°  ^ If  (k)  (t)  =0  for  k  =  0,1, . . . ,n  } 

»  |t|~ 

A  natural  norm  in  X°  is 

n 


a  t  l(k) 


f  II  =  l  sup  I  e”  1  ” 1  f 
k=0  tefi 


(t) 


Moreover,  let  Y°  =  {f  c  cn(  1R,  IR)  |  lim  e°  ^f  ^  (t)  =  0  for  k  =  1, . . .  ,n, 

n  Itl- 


lim  e  0tf  (t)  =  0} 
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A  natural  norm  in  Y° 

n 


is 


II  f  II  *  l  sup  |eCT^f  (k)  (t)  I  +  sup  |e-atf(t)|  +  suple^ffft)  -  f(»)| 
k=l  t  fIR  t<0  t>0 

+  |f {•)  | 

Theorem  2.4  Let  again  Y  denote  (y»  Y-.  *Y->»  •  •  •  »Y  ,6, ,...6  )  and 
XL  K  1  2  n  1  n 

VQ  =  ^et  0  >  0  be  small  enough  (smaller  than  all  the  real  parts 

i  Ai 

of  the  non-zero  eigenvalues  of  the  linearization).  Then  the  following  holds: 
If  f  e  X®  has  sufficiently  small  norm,  then  (2.2)  lias  a  solution  Y 
satisfying  Y  -  YQ  e  Yn+l  x  Y  depends  smoothly  on  f. 

PROOF:  When  we  put  Y  -  Yq  =  Z,  equation  (2.2)  can  be  written  in  the  form 

G(Z,f)  =  0,  and  G  is  a  smooth  mapping  of  (Y°  ,  x  (x°  ,)2N)  x  X° 

n+1  n+1  n  into 

(Xn)2N*1‘  Moreover»  ifc  can  easily  be  concluded  from  lemma  2.1  that  the 
linearization  DZG  (0,0)  is  an  isomorphism.  The  implicit  function  theorem 
implies  the  result. 


3.  GLOBAL  BEHAVIOUR  OF  SOLUTIONS  FOR  LARGE  f 

THEOREM  3.1  Let  f :  ]R  -*■  3R  be  continuous  and  such  that  lim  e  0tf(t)  =  0 

(a>0  small  enough),  f (t)  *  0  for  t  >  t  .  For  every  such  f,  equation  (2.2) 

0  K. 

has  a  unique  solution  satisfying  lim  y(t)  =  1,  lim  y  =  lim  6^  =  -j-  . 

t-*  -»  t-*-  t-+-  — »  i 

This  solution  exists  globally  for  all  times  t  and  lim  y(t)  exists  and 

t-*-  +00 

is  strictly  greater  that  zero. 

PROOF;  If  t^  is  chosen  large  enough,  e  0tf(t)  becomes  small  on  (-»,  -t^) 

and  one  can  use  an  implicit  function  argument  analogous  to  theorem  2.4  to 

prove  the  existence  of  a  solution  on  (-“,  -t^).  This  solution  is  unique  in 

the  class  of  solutions  approaching  the  limiting  values  at  t  =  -“>  at  a 
at 

rate  of  e  .  However,  if  a  solution  tends  to  these  limits  at  all,  it  can 
be  seen  from  the  last  two  equations  of  (2.2)  and  the  implicit  function 
theorem  that  y^  and  <5^  tend  to  their  limiting  values  at  a  rate  of 
e0t.  The  first  equation  then  implies  that  y  approaches  its  limiting 
value  at  the  same  rate.  Hence  the  solution  is  actually  unique  in  the 
class  of  all  solutions  approaching  the  prescribed  limits  at  t  =  -°°  as 
claimed  in  the  theorem. 

We  now  continue  this  solution  to  the  right,  and  we  have  to  make  sure 

that  it  does  not  blow  up  in  a  finite  time.  For  that  purpose  it  is  more 

convenient  to  consider  (2.1)  rather  that  (2.2).  From  the  second  and 

third  equation  we  see  that  as  long  as  y  stays  positive,  and 

have  a  positive  lower  bound  for  all  finite  times,  which  is  independent  of 

3  ot 

y.  Hence,  if  y  becomes  too  large,  g^y  will  dominate  over  fy  and 

also  over  h.  (since  this  is  less  than  some  constant  times  max  y(x)). 

1  (-°°,t) 

Analogously,  if  y  becomes  too  small,  will  dominate  over  fya  and 

g^y3.  Hence  y  cannot  go  to  zero  or  infinity  in  finite  time,  whence  we 

find  global  existence. 
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Let  now  t  >  tQ. 


Then  f  =•  0,  and  using  (2.2)  again,  we  find 


N  (  a. a. 

Ji  _ 

2  K, 


l 

i=l 


+  u 


Mi  1 


K. 

1  „  l 

a.  +  •: —  8.  +  t — 

l  A.  l  A. 

l  l 


N 

•  l 

i=l 


,  z 

A.p  a. 

gi  ] 

N  -s 

l  (a.  -  8.) 

i=l  1  1 

2  K.  +  XiP 

Ki  ! 

a.  +  -i 

1  A. 

1 

8.  +  -5— 

1  A. 

1  1  i 

(3.1) 


K. 


K. 


Here  we  have  put  a.  =  y.  -  ,  8.  =  <5.  -  t—  .  As  we  know  that  y.  and 

i  iAjiiA.  i 


S.  stay  positive,  the  denominators  a.  +  ,  8.  +  y-  are  always  positive, 

1  i  i 

and  the  left  side  of  the  equation  (3.1)  is  thus  the  derivative  of  a  positive 
definite  function  that  decreases  along  trajectories.  As  an  immediate 
consequence  we  obtain  that  cu  and  8^  tend  to  0  exponentially  for 
t  -*■  °°.  One  easily  concludes  from  (2.2)  that  y  approaches  a  constant. 
COROLLARY  3.2  If  f  is  always  non-negative  and  not  identically  zero,  then 
y  (<*>)  >  y  (-«>) ,  if  f  is  always  non-positive  and  not  identically  zero,  then 
y(°°)  <  y  (-00) . 


PROOF;  Assume  f  >  0,  the  other  case  is  analogous.  It  is  immediate  from 

the  integral  equation  (1.1)  that  f  >  0  implies  y  >  1  for  all  t. 

Moreover,  if  f  /  0,  there  must  be  some  t*  such  that  y(t*)  >  1. 

Let  now  z (t)  =  min  y(t).  Then  (1.1)  implies  that 
t  e  [t*,t] 


U 

f-37-)  z(t)  >  min  (0,-  —  /  a  (t-s)  (z3  (t)-l)ds) 

ldtJ+  U  _oo 


t* 


j  a  (t-s)  (z3(t)-l) 


If  z(t)  -  1  is  sufficiently  small,  this  gives  an  inequality  having  the 
form 
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a 


z  (t)  >  -  Ce  >Ct  (z-1) 


It  follows  immediately  that  t  z(t)  >  1. 

We  conclude  this  chapter  with  an  argument  showing  that  theorem  3.1 

does  not  hold,  if  the  condition  that  f(t)  =  0  for  t  >  tQ  is  replaced 

by  exponential  decrease  of  f  and  a  ^  1  (in  case  a  =  1  the  previous 

N 

argument  still  goes  through,  the  only  difference  being  that  f (t)  £  (a.  -8.) 

i=l 

has  to  be  added  on  the  right  side  of  (3.1)).  We  restrict  ourselves  to  the 
case  N  =  1.  Now  (2.1)  reads: 

•  3  ,  .  a 

-py  =  gy  -  h  -  f  (t)y 

K 

g  ■  -Xg  +  — 
y 

h  *  -Xh  +  Ky 

We  now  solve  these  equations  for  t  >  0  by  the  following  ansatz: 

vt  -2vt  -Xt  ,  vt  -At, 

y  =  yQe  ,  g  =  gQe  +  g1e  ,  h  =  hQe  +  h^e 


_  _  (l-a)vt  3-a  <(3-a)v-X)t  ~a  (-av-X)t 

f  =  V  +  giYo  6  -  hiyo  e 


After  some  calculation  one  finds  that  this  satisfies  the  equations  if 


-  K 

g0  2,,  .  ,  '  0  X  +  v 

yQ  (^-2v) 


_  a-1  3vK  +  pv (X-2v)  (X+v) 

f  y  E  — - - 

00  (X-2v) (X+v) 

We  thus  find  solutions  where  f  goes  to  zero  exponentially,  but  y  -*■  00 
for  a  >  1  and  y  -*■  0  for  a  <  1. 

All  we  have  to  make  sure  is  that  by  appropriate  continuation  for  t  <  0 
we  can  match  the  conditions  at  t  =  For  this  purpose  continue  y  in 

an  arbitrary  way  to  the  left  such  that  y  is  smooth  and  approaches  1 


-i.  j  W  'll.  SOaiMVUWMMr  .  ■ 


exponentially  at  t  *  The  equations  for  g  and  h  respectively  then 

have  unique  solutions  approaching  for  t  -*■  -<®.  These  solutions  can  be 
matched  to  the  solutions  for  t  >  0  by  appropriate  choice  of  g1  and  h^ 
Finally  f  is  determined  by  the  first  equation. 


4.  THE  CASE  y  ■  0 


In  this  case  the  first  equation  of  (2.1)  becomes 
,  N  N 

y  •  I  <3*  -  l  h  -  f(t)y“  =  0. 
i=l  1  i=l 

PROPOSITION  4.1  For  any  g  >  0,  h  >  0  and  0  <  a  <  3  the  equation 
F (y)  =  gy3  -  h  -  fya  =  0  has  a  unique  solution  in  (0,“). 

PROOF;  We  have  F(0)  <  0,  lim  F(y)  >  0,  so  there  is  clearly  a  positive 

y-K» 

solution.  To  show  it  is  unique,  we  investigate  zeros  of  F'  (y) .  We  have 

Ff(y)  *  3gy2  -  afyU  1.  If  y  >  0  and  F'(y)  =  0,  we  find 
1  3  3 

F(y)  =  ^yF* (y)  +  y  (1  -  —  )g  -  h  <  0.  This  means  F  cannot  have  a  positive 
maximum,  whence  the  result. 

The  solution  y(g,h,f)  can  then  be  inserted  into  the  other  equations, 
yielding  a  system  of  2N.  equations. 

THEOREM  4.1  The  same  statement  as  in  theorem  3.1  holds  also  for  u  *  0. 
Also,  Corollary  3.2  still  holds. 

SKETCH  OF  THE  PROOF:  The  existence  of  a  solution  on  (-",-t^  and  global 
existence  in  time  are  proved  in  the  same  manner  sis  before,  and  we  do  not 
repeat  the  arguments.  If  f  =  0,  one  now  finds  from  (2.2) 

Y .  =  -A.y.  +  K.  +  2yX 

i  ii  i  ly 

6.  =  -A. 6.  +  K.  -  6.^ 
i  ii  l  iy 

This  leads  to 


N 


fi  Vi  .  ®i8i  1  ? 

r  A.  a. 

1  ■ 1-  1 

A.  8.  \ 

1  1 

[2  K,  +  K.  J 

W  1=1 

1  1 

1  Ki 

2  (a,  +  -r^ 

1  A. 

1 

1 

+  y  •  i  "  V 

y  i=l 


where  and  8^  are  defined  as  before. 

Since  J(a^  -  8^)  is  now  equal  to  zero,  we  still  find  that  and  8^ 

approach  0  exponentially,  whence  the  result. 

-11- 


For  the  corollary,  observe  that 


y (t)  •  3y2(t)  /  a(t-s)  — =  -  /  a'(t-s)  - 
-00  y  (s )  -oo  v 


v3(t) 
ty2  (s) 


y  (s) 


ds 


Using  this,  one  can  apply  an  argument  analogous  to  the  previous  one. 
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